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Abstract 

The self-similar equilibrium models of self-gravitating, rotating, isothermal systems are investigated 
analytically. In these models the rotation velocity is constant and the density varies as 5 , where r and 
6 are the spherical radius and the co-latitude, respectively. The nonaxisymmetric solutions contain three 
free parameters, one of the parameters depends on the rotation velocity. These parameters determine the 
overall shape of the density distribution. By assuming that the dominant component of the magnetic field 
is purely toroidal and the ratio of the purely toroidal magnetic pressure to the gas pressure, a, is spatially 
constant, the axisymmetric solutions generalized so as the effect of magnetic field could be studied. We find 
that the equilibria of axially symmetric systems yield ellipsoids or spheres only when the ratio of rotation 
velocity to the sound speed is taken to be \/2a. 

Subject headings: hydromagnetics - stars: kinematics - stars: formation - galaxies: structure 



1. INTRODUCTION 

Understanding the equilibrium structure of self-gravitating isothermal systems is a fundamental issue in astro- 
physics. In star-forming regions many structures have been identified from sheets, to filaments, to elongated 
clumps (Kulkarni & Heiles 1988; Jijina, Myers, & Adams 1999). By some analytical solutions, the hydrostatic 
equilibrium of these structures have been studied previously (Shu 1992; Ostriker 1964). Several authors, how- 
ever, investigated the problem by numerical techniques so as to include the effects of rotation and magnetic 
fields (see, e.g., Fiege & Pudritz 2000; Curry 2000; Tomisaka, Ikeuchi, & Nakamura 1988). Many models of 
newborn stars built upon these equilibrium configurations (e.g., Kiguchi et al. 1987). Also, there is a close 
analogy between gaseous self-gravitating equilibrium objects and collisionless stellar systems (see, e.g., Binney 
& Tremaine 1987). Thus the equilibria of isothermal systems can be used for determining the steady-state 
distribution function of collisionless stellar systems (Toomre 1982; Hayashi, Narita, & Miyama 1982). 

With these considerations, recently Medvedev and Narayan (2000; hereafter MN) classified and derived 
analytically self-similar axisymmetric equilibria of a self-gravitating isothermal system. They presented two 
families of equilibria: (1) Cylindrically symmetric solutions, in which all quantities depend on the cylindrical 
radius only. (2) Axially symmetric solutions, in which the quantities vary as functions of the spherical radius 
and the co-latitude. In fact, these axially symmetric solutions are very interesting. These equilibria form a 
two-parameter family of solutions. One of the parameters depends on the rotation velocity and the other one 
controls the symmetry of the solution with respect to the equatorial plane. Nevertheless, this study like most 
previous work has focused on the case of axisymmetric equilibria of self-gravitating systems. The absence in the 
literature of nonaxisymmetric equilibrium solution suggests that a fresh theoretical approach is necessary. 
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Our motivation for this work is to determine nonaxisymmetric equilibria of self-gravitating systems. We 
extend the work of MN by obtaining nonaxisymmetric self-similar equilibria of a self-gravitating system ana- 
lytically. Because of great complexity and intrinsically nonlinear nature of the problem, any analytical solution 
is unique and, thus, of great scientific value. Our nonaxisymmetric equilibrium configurations form a three- 
parameter family of solutions and in the case of axisymmetric equilibrium, these solutions reduce to the solu- 
tions of MN. Furthermore, we shall study the effect of magnetic fields on the axisymmetric self-similar equilibria 
solutions. Since we are interested in treating the problem analytically and for simplicity, we assume that the 
dominant component of the magnetic field is purely toroidal and the ratio of the magnetic pressure to the gas 
pressure, a, is spatially constant. MN showed that the non-rotating axially symmetric solutions form confocal 
ellipsoids or spheres. But we find that the non-rotating, magnetized, axially symmetric solutions tend to the 
cylindrically symmetric structures. However, the equilibria of axially symmetric systems both with rotation and 
toroidal magnetic field yield ellipsoids or spheres only when the ratio of rotation velocity to the sound speed is 
taken to be V2a. 

In §2, we present the general formulation and basic assumptions. We derive a three-parameter family of 
nonaxisymmetric self-similar solutions in §3 and investigate the properties of these solutions in §4. We summa- 
rize our results in §5. 

2. GENERAL FORMULATION 

We present here the basic equations used to describe the equilibrium structure of a rotating self-gravitating 
system permeated by a toroidal magnetic field. We assume that the medium is an isothermal ideal gas and 
the magnetic field is frozen into the gas. However, the assumption that the system is isothermal is good for 
low-density region. In spherical coordinates (r,9,ip), the basic equations are magnetohydrostatic equation, the 
Poisson equation and the equation of state: 
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p = c 2 s p, (5) 

where p, p, v v , c s , B v and denote the gas density, pressure, toroidal component of the velocity due to the 
rotation, sound speed, toroidal magnetic field, and gravitational potential, respectively. We assumed that the 
entire system is rotating around a common axis which specifies the axis of the spherical coordinate system. 
It is clear that if Q, denotes the angular velocity, then v v = ilrsinO. We note that purely toroidal magnetic 
field under the cylindrical symmetry, B v (r, 9), automatically satisfies V.B = 0. In fact this configuration of the 
magnetic field requires a rather artificial current configuration. Namely the current must flow along the rotation 
axis (z-axis). Although systems with such fields are not too common in astrophysics, our results provide an 
intuitive way to understand more complex systems. Also, the rotation velocity, the density and gravitational 
potential, in general, are functions of r, 9, and ip. 

Solving equations (1) to (5), in general, is a difficult work. Thus to simplify the problem, we need another 
constraint. As for molecular clouds, observations suggest that the magnetic field strength often varies with the 
density roughly according to B ~ p^ (e.g., Heiles et al. 1993). It implies that the ratio of magnetic pressure to 
the thermal pressure, a, is spatially constant: 

where in molecular clouds observations show that < a < 10 (Heiles et al. 1993). For our problem we use it as a 
free parameter. By changing the value of a, the effect of purely toroidal magnetic field on the equilibrium state of 
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self-gravitating systems can be studied. It must be noted that by using this constraint we can investigate effects 
of magnetic field only on the axisymmetric solutions. In fact our nonaxisymmetric solutions are nonmagnctized. 
Now, we introduce the dimensionless variables according to 

p -> pp,r -> rr, -> ^ -> vv v , B v -> BB V , (7) 

where 

p = p ,r = (47rGp yic 2 s ,i> = c 2 s ,v = c Sl B = (4Trp Q )i Cs . (8) 

By transforming to dimensionless variables and using equation (6), equations (1) to (4) are rewritten as the 
following: 
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By substituting equations (9), (10), and (11) into equation (12) we obtain 

a + 1 d 2 dlnp a + 1 d Q dlnp 1 d 2 lnp 

It would be difficult to study systematically the very large space of all possible solutions of these equations. In 
fact, we need some information about the density, p, and the rotation velocity, v v . For example, by knowing 
the form of one of the variables it is possible to find another variable. But here we are interested in self-similar 
solutions which will be studied in the next section. These solutions are interpreted with much less effort than 
would be involved in a direct attempt to solve the full set of partial differential equations. 

3. SELF- SIMILAR SOLUTIONS 

We introduce the self-similar forms for the density and the rotation velocity as 

p(r,d,ip) = )«(p = 7 — r-rr; — , (15) 

v ; r"i v (rsin6) v * 

where v 2 , and f(0,<p) to be calculated. The same self-similar solutions for the axisymmetric case has been 
introduced by MN. Their angular part of the density distribution was a function of 9 only. Since we are going to 
study the nonaxisymmetric solutions, the angular part of the density is a function of both 9 and <p. Substituting 
the self-similar forms into equations (13) and (14), we obtain 

v 1 =2,v 2 = 0, (16) 

and f(9, <p) can be found by solving this equation: 

- 2{a + + ^9d9 ismd ^9- ) + ^9^ +f = °- (17) 

This is the main equation which must be solved. We note that for the nonmagnetized and axisymmetric 
case, this equation reduces to the equation which has been solved by MN. To solve equation (17), we use 
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transformations which are equivalent to Howard's transformation used by Toomre (1982) and MN. First, we 
rewrite this equation as follows, 



19. . .9, ,, . 2 ... 1 dHnf 



[sinO— ln{f sin 2 0)] + -_ — -f + f = 0. (18) 
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Upon introducing a new function, w, and a new independent variable, £, as follows, 



9 1 

w = ln(fsin 2 6),£ = ln\tan{-)\^+' 1 , (19) 



equation (18) is greatly simplified and becomes 

d 2 w d 2 w 

w + w +eW=0 - m 

This equation is similar to the Lanc-Emden equation which has been written in cartesian coordinates. 
Schmid-Burgk (1967) and independently Stuart (1967) presented an interesting two-dimensional solution of the 
Lane-Emden equation in cartesian coordinates. Thus, in analogy to their solution we can write an analytical 
solution for equation(20), 

-7 r 2q 2 (1-C 2 ) 1 

W ~ 1 [cosh(a£ + b) + Ccos(ap>)} 2 } ' [ ' 

where a, b, and C are free parameters. Recalling the transformations (19), we finally obtain an analytical 
solution, 

lieM i ^2 . m 

sm [cosh(ln\tan^\V^ + b) + Ccos{ap)] 2 
This solution can be simplified further and the densiy becomes 
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where A = 2a, B = e 2b , and — 1 < C < 1. This solution describes the nonaxisymmetric density distribution of a 
self-gravitating system. It is clear that to obtain the nonmagnetized axisymmetric equilibria we must set a = 
and C = and our solution reduces to the solution of MN. Equation (23) presents a three-parameter family of 
solutions (note that for nonaxisymmetric we must set a = 0) and is smooth and well-behaved for all values of 9 
and ip in the domains < 9 < tt and < ip < 2n. However, the equation is ill-defined at r = or 9 = 0. We can 
investigate if a singular solution exists at the origin or on the axis. Using Gauss'theorem, MN calculated these 
additional mass densities for the case of the nonmagnetized axisymmetric equilibria. The reader is referred to 
MN for a full derivation. Following MN, we find that our solution has no singular mass at the origin. But the 
requirement of a non-singular density on the axis gives, 

\/a + 1 

This equation shows that the rotation velocity and the magnetic field determine the value of A. Thus, the 
solution (23) becomes 
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Since the exponent of tan| determines the overall shape of the density distribution, equation (26) is very im- 
portant. In the next section we study the properties of the solutions. 

4. PROPERTIES OF THE SOLUTIONS 

Since a real self-gravitating system is finite, we truncated the system by an external pressure. The external 
medium is assumed to be non-self-gravitating and of negligible density. In this way we can define the surface 
of the system. If p s represent the external pressure in a nondimentional form, we obtain the equation of the 

surface 

2A*B(1-C*) 1 \tan\\^f (2?) 

Ps sinO 1 + B \ tan e_\A cS +2CVB\tan^co S ^ 

Notice also that for a — and C = 0, this equation reduces to the solution which MN obtained. As equations 
(24) and (26) show, the rotation velocity, vq v , and the ratio of magnetic pressure to the gas pressure, a, determine 
the values of A and A e f. But it is possible to make another constraint: Since we are interested in the solutions 
which have no discontinuity on the surface, we accept those values of A that r s (0, tp = 0) = r s {9 1 ip = 2n). Thus, 

A = 2u, (28) 

where v is a positive integer. So only a limited number of values of v 0lfi and a gives acceptable values of A. 

As mentioned before, the value of A e f determines the overal shape of the system. For the nonmagnetized 
case, we have A e f = A. So, only the rotation velocity determines the value of A. As the v 0ip increases, the 
value of A becomes larger. Equation(26) shows that as a increases, the value of A e f decreases. It seems that a 
purely toroidal magnetic field with a constant ratio of the magnetic pressure to the thermal pressure essentially 
cancels the effects of rotation. With these considerations, we focus this current discussion on nonmagnetized 
and magnetized solutions. 



4-1. Nonmagnetized Equilibria 

First, we study the nonrotating solutions: A e f = A = 2. The solution (27) reads 



' ^' ^ 1 + ecos9 + 8sin9cosip ' 



where 



8(1 -CP) 1 - B 2CVB 

= (STf' e = TT~b T+W' (30) 

Note that the case C = and B ^ 1 corresponds to a prolate ellipsoid and if B = 1 this ellipsoid reduces to 
a singular sphere. These axisymmetric equilibria have been studied extensively by MN. They showed that the 
solutions with B > 1 determine prolate configurations which are shifted upwards along the z axis and those 
with B < 1 arc shifted downwards. So, the parameter B controls the symmetry of the solutions with respect to 
the equatorial plane. 

It can easily be verified that the general case C ^ and B ^ 1 represents a family of ellipsoids. By rotating 
the coordinate system (x, y, z) about y axis from the positive direction of z axis to the positive direction of x axis 
we can define a new coordinate system (X, y, Z). Assuming that ui is the angle of rotation, upon straightforward 
but cumbersome algebric manipulations, we can rewrite equation (29) in a new coordinate system (X,y,Z) 

where ^ ^ 

— esinuo — Scoslj, 7 = ecosu; + Ssinuj, u) = -tan^ 1 -= -r, (32) 
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Equation (31) shows that the case C = and B ^ 1 represents a family of ellipsoids, with traces on the 
y — z and x — z planes which are ellipses and traces on the x — y plane which are circles with centers at origin. 
Since r s {9) is not invariant under the transformation 9 — > 7r — 9, this configuration is not symmetric with respect 
to the equatorial plane. On the other hand, the case B = 1 and C / corresponds to a family of ellipsoids 
which are elongated along the x axis and r s (9, ip) is invariant under the transformations 9 — > ir — 9 and ip —* ip. 
Thus, the equatorial plane is the symmetry plane and projections of these ellipsoids onto the y — z plane are 
circles with centers at origin. In fact, the general solution (27) shows that r a (6, ip) is invariant under these 
transformations only for B = 1. So, the parameter B in our nonaxisymmctric solution controls the symmetry 
of the solutions with respect to the equatorial plane. 

As mentioned before, we restrict our study only to the solutions which have no singular mass on the axis. 
This constraint gives A y 2 or v > 1; note that the case v = 1 has already been studied. For C = and A > 2, 
equation (27) corresponds to a family of toroidal structures in which the parameter B determines the asymme- 
try of these equilibria with respect to the equatorial plane. In this case, as A increases, the systems become 
flattend and tend to a disk. It seems that in our solution the parameter C alone determines the asymmetry of 
the system with respect to the z axis. Although it is true, the parameter A (or v) is also very important. We 
see that the equation (27) is invariant under the transformations 9^9 and ip — > n + ip only for even values of 
v. Therefore, in the case C / we may have configurations which are symmetric with respect to the z axis. 
The rotation velocity has a fundamental role in determining the overal shape of the density distribution. For 
small values of C and A > 2, the toroidal equilibrium structures are elongated in some directions. We refer 
to these elongated parts to as knobs. As discussed above, the number of these knobs depends on the rotation 
velocity. For example, the case A = 8 [y = 4) represent a toroidal configuration with four knobs. With increas- 
ing C and tending the parameter to the unity, these knobs become larger and the toroidal part becomes smaller. 

4-2. Magnetized Equilibria 

It is useful to study the effect of the magnetic field on our solutions. We postulated that the magnetic field 
structure is purely toroidal and corresponds to that of constant ratio of the magnetic pressure to the thermal 
pressure. Although purely toroidal fields are probably far from reality, this is the best way to investigate the 
effects of the magnetic fields on the equilibria analytically. 

Now, we consider the effects of the magnetic field on axisymmetric solutions: C = 0. Equation (26) 
shows that a nonrotating system always has A e f < 2 and the value of A e f decreases with increasing the 
strength of magnetic field. While the nonrotating and nonmagnetized systems are ellipsoids, the nonrotating 
and magnetized systems tend to the cylindrically symmetric configurations. Also, if we consider rotation such 
that Vo v < V^ce, the system has cylindrically symmetric structure. 

If voip = V2a, we have A e f = 2 and the equilibria are like those which have been investigated in the previous 
subsection for the case C = 0. Finally, the case vo v > V2a and C = corresponds to a family of toroidal 
configurations. We know that as vo v — > oo the equilibria tend to a disk. However, the magnetic field decreases 
the effect of rotation and in this case it causes the system to become a toroidal structure not a disk. 

5. DISCUSSION AND SUMMARY 

Self-similar equilibria of a self-gravitating rotating system containing a purely toroidal magnetic field has been 
investigated. A three-parameter family of solutions for nonmagnetized self-gravitating systems have been found. 
By assuming that the ratio of the magnetic pressure to the gas pressure, a, is spatially constant, the axisym- 
metric solutions were generalized so that the effect of the magnetic field could be studied. We have shown that, 
depending on the values of these parameters, the overall behavior of the density distribution changes. 

Recently, Galli et al. (2001) studied binary and multiple star formation by considering the nonaxisymmetric 
equilibria of self-gravitating, magnetized, differentially-rotating, completely flattened singular isothermal disks 
with critical or supercritical ratios of mass-to-flux. They found that lopsided configurations exist at any di- 
mensionless rotation rate. Also, multiple-lobed (M — 2,3,4, ...) configurations correspond to rotations of the 
equilibrium configurations by multiples of jj and bifurcations into sequences with M — 2,3,4,5, and higher 
symmetry require considerable larger rotation rates. These results are in a well agreement with our analysis, 
at least qualitatively. We must note that the flattening of our nonmagnetized equilibria produced by rotation 
rather than by magnetic fields. 
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As equation (25) shows the parameters A or A e f determine the overall shape of density distribution. We 
showed that only a limited number of values of vo v and a give acceptable values for A or A e f. Since for 
nonmagnctized solutions we have A — A e f, as the value of rotation velocity increases, the value of A becomes 
larger. Equation (25) shows that for A > 2, the toroidal equilibrium structures have knobs and the number of 
these knobs depends on the rotation velocity: as rotation velocity increases, the number of the knobs increases. 
It is interesting that in the case C / 0, some values of A construct equilibrium configurations which are 
symmetric with respect to the z axis. We sec that our multiple-lobcd (y = 1,2,3, ...) equilibria require larger 
rotation rates. The equilibria solutions become flattened as the rotation velocity increases and they tend to 
a thin disk as A — > oo. Thus, our nonaxisymmetric solution tend to a nonaxisymmetric thin disk as the 
rotation velocity increases. These configurations are equivalent to the nonaxisymmetric solutions which have 
been obtained by Galli et al. (2001). Detailed observations of the gas disks associated with bipolar outflows 
indicate that the disks arc rotating very fast (Kaifu 1987). The rotation period is comparable to the free-fall 
time scale of the gas disks and it means that the disk is near the rotational balance. Thus, in these systems the 
effect of rotation should have been important in the formation of the central star. On the other hand, it is very 
important that the formation of multiple stellar systems could never result from any calculation that imposes 
a priori an assumption of axial symmetry. So, as long as the starting conditions are the nonaxisymmetric disk 
solutions, gravitational collapses seem in general to produce fragmentation. 

For B — 1 and C = 0, the solutions have been obtained by Toomre (1982) and Hayashi, Narita, & Miyama 
(1982). In fact, Toomrc-Hayashi model is a good approximation to a centrifugally supported disk formed by 
collapse of a rotating gas cloud (see, e.g., Nomura & Mineshige 2000). Since the evolution of such systems 
may be driven by instabilities, Hanawa, Saigo, & Matsumoto (2000) discussed the stability of the Toomre- 
Hayashi model against axisymmetric perturbations and Shu et al. (2000) analyzed the stability of this model 
under thin disk approximation while taking account of magnetic fields. The present analysis suggests that 
it is necessary to study the stability of our nonaxisymmetric disk solutions. Of course, the nonaxisymmetric 
solutions are nonmagnetized and the self-gravity balances with centrifugal force and gas pressure. But the 
magnetic field reduces the effective gravity and increases the sound speed (see, e.g., Li & Shu 1997). The 
effects can be fully taken into account by replacing the gravitational constant (G) and sound speed (c s ) with 
effective ones. So, the stability of magnetized disk is the same as that of nonmagnetized one. Also, one should 
perturb the self-similar equilibria and investigate their stability. The possibility of gravitational instability and 
its nonlinear consequence should be considered. Particularly, the nonaxisymmetric equilibria may fragment due 
to the gravitational instability and multiple stellar systems may form. 

Since magnetic field, as well as rotation, plays an important role in star formation, an effect of purely toroidal 
magnetic field on the axisymmetric solutions has been studied. We assumed that the ratio of magnetic pressure 
to the thermal pressure, a, is spatially constant, since we found this case to be integrable analytically. We 
showed that nonrotating and magnetized systems tend to the cylindrically symmetric configurations. On the 
other hand, MN showed that the nonrotating and nonmagnetized systems form ellipsoids (or spheres). If a 
system rotates such that vq v < y/2a, the system has a cylindrically symmetric structure. But if vq v — ^/2a, 
the systems form ellipsoids or spheres. For vo v > V2a the axisymmetric equilibrium corresponds to a family of 
toroidal configurations. In fact, the magnetic field decreases the effect of rotation. 

MN proved that the axisymmetric solutions with B ^ 1 are not force-free at two singular points, r = and 
r = oo. In fact external forces hold the axisymmetric system in equilibrium. Are the nonaxisymmetric equilibria 
likely to be not force-free at singular points? Although we did not investigate this problem, it is likely that 
the nonaxisymmetric solutions are also not force-free at singular points. However, we think such configurations 
represent legitimate states of equilibrium. In this regard, we would like to remember the explanation of Galli et 
al. (2001): "in the solutions, the (infinite) gravitational force at the origin is exactly balanced by an (infinite) 
pressure gradient acting in concert perhaps with an (infinite) centrifugal force. This balance is qualitatively no 
different than at any other point in the system , and it would be an artificial restriction to rule out eccentric 
equilibria simply because they have a nontrivial balance of forces at the origin rather than a trivial one." 
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